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Abstract 

An expression, suitable for practical calculations, for the Green function G{E; xi, X2) of 
the Schrodinger equation is obtained through the approximations of the path integral by n- 
fold multiple integrals. The approximations to ^{G{E;x,x)} on the real £^-axis have peaks 
near the values of the energy levels Ej, where exact di{G{E; x, x)} behaves like S{E — Ej). 
The analytic and numerical examples for one-dimensional and multi-dimensional harmonic 
and anharmonic oscillators, and potential wells, show that median values of these peaks for 
approximate G{E; 0, 0) corresponds with accuracy of order 10% to the exact values of even 
levels already in the lowest orders of approximation n = 1 and n = 2, i.e. when the path 
integral is replaced by a line or double integral. The weights of the peaks approximate the 
values of the squared modulus of the wave functions at a; = with the same accuracy. The 
accuracy depends on the type of the paths used. 



* Preprint ITEP-111-82 



1 Introduction 



Continual (functional) integrals, or path integrals, introduced to quantum mechanics by R.Feyn- 
man ^ are widely used for constructing series in perturbation theory, for finding quasi- 
classical asymptotes, for quantizing gauge fields [Sj-IZI- Mathematical questions of the functional 
integral theory are discussed in Path integrals combined with Monte-Carlo technique are 

being used for solving problems of quark confinement jlUj-|12j. 

However, the question of applying path integral approximations to problems of practical 
computations of quantum system spectra remains virtually unexplored. First, not very success- 
ful, attempts of this kind are described by Kac|13j. An example of computation for a system 
of two particles sitting in a potential well is considered in Among more recent papers it 
is worth noting the work ^S], where an estimate for the ground level from below is obtained 
with a help of path integral. Papers ^Bl-lttSl employ "Wick rotation" and they allow obtaining 
information only about the ground level, moreover, for some examples considered in the 
results disagree with exact values by a factor more than two. It seems that path integrals were 
applied for global spectrum analysis in quasi-classical limit only ^HJ-CHl- It remains unclear if 
approximations to path integrals can have a practical value for calculations of the Schrodinger 
equation spectra. To clarify the situation, it seems appropriate to begin with simplest quantum 
systems and to investigate the influence of the path types and of the order of approximations 
on the results. 

Here we suggest the following approach to finding the energy levels Ej of a quantum sys- 
tem with the help of path integrals. The kernel of evolution operator IC{xb,T; Xa,0) (i.e. the 
time-dependent Green function) can be expressed via eigenfunctions ^jix) of the Schrodinger 
equation: 

IC{xb,T;Xa,0) =y2'4'j{xb)ipUxa)ex.p{-iEjT/h). (1) 



j 



For E shifted for an infinitesimal step to the upper half of the complex plane and for Xa = Xh 
xq, equation (1) implies: 



n 



(2) 



where G{E,xo) = G{E, xq,xo) is the Green function of a stationary Schrodinger equation. 

On the other hand, the kernel (1) can be expressed via the path integral in configuration 
space dH]: 

IC{xb,T;Xa,0) = lim /C„(xb,r;xa,0); (3) 



n— >oo 
-oo 



/C„(x6,r;xa,0) =C„F(fr)-("+i)/2 J l[dxjexp[ii/h) J L{qF,qF,t)dT\. (4) 

-oo 

Here 

qF{t)=q*{t)+XF{t), (5) 

with (7*(t) - a segment of a straight line connecting Xa and x;,, while XF{t) is a broken line 
with vertices in Xj {j = 1, 2, ... n) separated by equal time intervals T/ (n + 1); L{qF, Qf, t) is a 
Lagrangian, and 

C„i. = [m(n + l)/(27rn)]("+i)/2^ (6) 
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Subscript F in qp and CnP denotes the type of the path xp{t) (Feynman). 

Let us define a function Gn{E,xo), an approximation of the n-th order to the exact Green 
function G{E,xq): 

oo 

Gn{E,xo) = J )Cn{xo,T;xo,0) exp{iET/h)dT. (7) 



For n — oo the real part of Gn{E,xo), i.e. 3fi{G„(£', xq)}, tends to a sum of ^-functions in 
the r.h.s. of (2), according to (3). For finite n (for a sufficiently good approximation) the real 
part ?R-{GniE , xq)} must have maxima near exact values of Ej. 

We put forward a technique for computations of Gn{E, xq) and carry out them for harmonic 
and anharmonic oscillators in one- and multi-dimensional cases and in a one-dimensional poten- 
tial well of a shape Uq cosh~^(g/a). We find a number of analytical expressions for 3?{Gi(-E', xq)}. 

An unexpected result of our numerical experiments is the fact that a "good" approximation 
for finding the values of energy levels (with accuracy of order 10 per cent) is Gn{E,xo) already 
for n = 1 and n = 2. 



2 Method of calculation 

The approximations /C^ in (4) can be built based on an arbitrary complete set of functions, not 
necessarily broken lines. For example, using an expansion over an orthogonal set of sines one 
can replace the paths qrit) in (4) by 



qs{t) = q*{t) + ^Xj sm{TTjt/T), 

3=1 



and the factor CnP by 



Cn 



n! 7r"/2 ^my"+i)/2 



(5.) 



(27r)V2 2" \nj 

The factor Cns is determined from the condition that an approximate expression for the kernel 
coincide /C„ with the exact /C for all values of n in case of a free particle. 
For a Lagrangian 

L{q,q,t)=mf/2-Uov{q/a), (8) 
with Uq and a being units of energy and length, we get from (4) and (5): 



oo 

GniiE,xo) = a''Cni j dTiiT)' 



(n+l)/2. 



W dx j exp 



-UoTfe{xo,Xi,...,Xn)+ET 



(9) 



Here the subscripts i = F, s correspond to paths qrit) and qs{t) respectively. All distances Xj 
and qi in the Eq.(9) and herafter are measured in units a. The following notation is used: 



j=0 



Xn+1 = Xo, 



niacin + 1) 



(10) 



2 



CTns = Y.j^^]^ = 7r2maV4; (11) 

i=i 

1 

fi(xo,xi,...,Xn) = j f{qeiT))dr, T = t/T. (12) 



Integral over T in (9) can be done analytically interchanging the order of integration over T 
and Xj. For E shifted to the upper semi-plane we have, according to [19], 

- oo i=i 

Here B^i = ira^Cne; H(n-i)/2(^) Hankel function of the first kind; 

z = li{onl[E' - Ft{xQ,xi, ...,x„)]}^/2^ (14) 

where E' = E/Uo, li = i2/h)iPiUo)^/\ 

The expression (13) is valid also as an approximation of the n-th order of a general Green 
function G{E;Xa,Xb). Then xq = Xa, Xn+i = Xf, in (10), and functions fi in (12) and (14) 
depend both on Xa and x^,. 

A useful tool for searching the levels can also be a spectral function, 

+0O 

Y{T) = I K:{x,T;x,0)dx = Y^gjexp (-^) , (15) 

— oo 3 

where gj is the degeneracy of the j-ih level. Defining 

oo 

F{E) = J Y{T) ex.p{iET/h)dT, (16) 



we see that Fne{E), i.e. the n-th order approximation to F{E) over paths of type i is given by 
the expression 

+°° n 

Fne{E) = aB„e f ncix,-(zK^)("-i)/2i7«_^/2(^)- (17) 
- oo ^=0 

This expression differs from (13) only by raising the order of integration by unity. Maxima of 
the real part of Fni{E) must lie near the exact energy levels while the height of the maxima 
must reflect the degree of the level degeneracy. 

The real part of the integral in (13) can be written in the form 

/n 
n tia;j(zK^)("-i)/V(„_i)/2(2:), (18) 

D 

where J(n^i)/2(^) is Bessel function. The domain D of integration over Xj in (18) is determined 
by the condition 

D : fe{xo,xi,...Xn) <E'. (19) 
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If the inequality (19) is violated then the integrand in (18) is purely imaginary. 

For even n = 2k + 2 (fe = 0, 1, 2 . . .) the Bessel function J(n-i)/2(^) reduces to elementary 
functions and we have from (18): 



D 



(20) 



The expression for real part of Fni{E) is similar to (18) and (20) with 11^=1 replaced by 

The formulae (9 - 18) can be easily generalized for a multi-dimensional case. If the dimension 
of vectors x and qis u then 



oo 



J n '^^j exp{(z/a)[/?£C7„^^//r - UoTfi^{xo, xi, f„) + ET]} 



_-oo J=l 



i=o 



cr 



E-2 ->2. 
J Xj ) 

j=0 



'e^{xo,xi,...Xn) = / ipiqeiT))dr 



(9a) 
(10a) 
(11a) 

(12a) 



+00 



r ^ J- i/(n + l) , 

Gn^i.(£^ ) = Bniu / 11 C?Xj(2;^/(T„^^) 2 X , (-^i.); 

B„,, = 7ra'^-C^,(2/?,/n)i-^; 

/n 
-| — r /^(ri+1) 1 

[[dxj{z^/aneu)~ x J ^(n+i) ^jz^). 

D ^=1 



(14a) 

(18a) 



In this work the calculations were done using formulae (18), (18a) and (20). We omit the 
subscript v below when i/ = 1. 



3 Results of calculations 
3.1 Harmonic oscillator 

The exact expression for the kernel fC of one-dimensional harmonic oscillator is well-known, still 
it is interesting to consider the behavior of Rne{E) for low orders of approximation. For 

'P{q) = g^ Uo = ni^/2'K (21) 
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it follows from (18) that, if n = 1, 



RiiiE') = Aee/'Ji/4i^e)J-i/4m, (22) 

where 

As = V2^, 6 = E'/4. 

In case of evenly symmetric potentials, U{q) = U{—q), for odd levels, the wavefunctions 
vanish at g = 0, ipji^) = 0; therefore, according to (2) maxima of the r.h.s. of expression (22) 

must correspond only to even energy levels Ej for j = 0, 2, 4, . . .. Those exact values are shown 
by arrows in Fig.l where the dependence of Rif and Rig on E' is plotted according to (22). 



1.5 




Figure 1: Dependencies Rii{E') for the one-dimensional harmonic oscillator (21). Solid £ = s, 
dashed £ = F. The exact energy values of even levels are shown by arrows. 

One can see that the dependence Ru{E') has clear peaks located in correspondence with 
the levels of the oscillator. E.g., for expansion over sines (5s) the location of the first maximum 
deviates from the exact value of the ground level by less than 15%. 

The behavior of Rni{E') (which is proportional to ?R.{Gni{E')} ) must reflect the properties 
of a sum of (5-functions. The location of a (5-function is determined by an integral relation. 
Therefore it is more natural to expect that a level corresponds not to a maximum of Rne{E'), 
but to an average of the integral (i.e. a median value) of Rni{E') between two minima of the 
dependence Rnt{E'). We call this median value a median of the peak, and the area under the 
curve 3f?{Gn£ (£■')} between two adjacent minima of Rni{E') will be called a weight of the peak. 
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The median of the ground level peak calculated from (22) is equal to 3.08 which differs form 
the exact value, tt, by 2%. 

For S> 1 we have from (22): 

Ru{E') « A^7r-iC7'/'[2-V2 + cos(2e^ - tt/2)], (23) 

i.e. maxima and medians of the peaks Rip(E') arc at Ej = (tt"^ + 1/2)), and of Ris{E') 
at Ej = 27r(j + 1/2), for large even j. Thus, medians of peaks of Ris{E') coincide with exact 
values of oscillator levels for E' S> 1, while the expansion over broken lines gives a spectrum 
which is qualitatively correct, but with a relative error ~ 10% (approximately the same as for 
the ground level). 

With account of the coefficient in (13), we obtain from (22) for i = s that the weight 
of the peak of the ground level is equal to 1.776{hmujy^^ . This differs from the exact value, 
7r^|V'o(0)P = {■Khmu)^/'^, by less than 1%. 

In general case: 

l^^'^^')!' = 0F22^(fc!)2 n ' 

therefore, we find an asymptotic formula 7r^|V'2fc(0)p = {irhmui/ky^'^ for levels with large number 
j = 2k. According to (23), for ^ = s the weights of peaks coincide with this exact value. For 
£ = F the difference between the peak weight and the exact value vr?i|-(/'2fc(0)P is of order ^ 10%. 

Now let us turn to approximations of the second order (n = 2). The results of numerical 
calculations of R2f{E') and R2s{E') based on the formula (20) are shown in Fig.2. For n = 1 the 
graphs of Rif{E') and Ris{E') were similar to each other, but now we see a drastic difference in 
the behavior of R2e{E'). The function R2f{E'), built with the use of broken paths, has "false" 
maxima, while R2s{E') demonstrates quite a regular behavior showing a tendency to a sum of 
(^-functions. 



2 (2^)! "iw 




5 10 5 10 

E'/2tt E'/2tt 

Figure 2: Dependencies R2e{E') for the harmonic oscillator potential (21). Left: for £ = F; 
right: for ^ = s. The exact energy values of even levels are shown by arrows. 
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It is likely that good results obtained with the paths (5s) for n = 2 are due to the fact that 
they are optimum for approximations of path integrals over the Wiener measure [20, 21]. For 
other potentials, e.g. for a box with reflecting walls, the behavior of the function RnpiE') is 
quite regular (see below) and it has no "false" maxima. Hence, the system of optimum functions 
should depend on the shape of the potential well. 

3.2 Anhcirmonic oscillators 

For power-law potentials 

V{q) = q"" (24) 
and even N we present the results of computations of Rnt{E') for n = 1, 2. Instead of E', here 
it is convenient to introduce a variable 

y = £;'(^+2)/2^(7r2ma2c/o/?i2)i/2/^i/^, (25) 

where for = 2p : p = 1, 2, 3 . . . the expression for Ijv reads In = (2p — l)!!/(2p)!!. The 
graphs R2s{y) for iV = 4, 10, 50 are shown in Fig.3. 

The case N = oo (i.e. a rectangular well of infinite depth) allows for an analytical expression 
for Ru{E). It follows from (18) that 

Ru{E) = 2777^ / Mx)dx, (26) 



where rip = 2\/2y/7r, rjs = y = 7r(i?m)^/^a/?i, and 2a is the width of the well. 

The integral (26) can be expressed through the Struve functions, yet it is more convenient 
to use tables and expansions for this integral given, e.g., in [23]. 

For rji ^ 1 , that is for E ^ Eq, expression (26) implies 

Ru{E) « 277^1 [1 + (2/7r?7^)V2 cos(% - 37r/4)]. (27) 

Maxima (and medians) of peaks of Rif{E) are reached at £^2fc^ = h{2k + 3 / 4) / VSma^ , 
whereas the exact values of energy levels number 2k are eH^ = h{2k + 1)/V8 ma?, for k = 
0, 1, 2.... The results of numerical calculations for the peak medians yii? and yis, corresponding 
to the ground level for oscillators with various values of N, are given in Table 1 for n = 1, 2. 

The data in Table 1 show that the difference between approximate and exact values for the 
ground level increases with growing N for d. = s and decreases for i = F. 

Medians of the peaks for higher levels also have a systematic shift relative to the exact values. 
The relative error for high levels grows from 0% for TV = 2 up to about 10% for A/" — oo when 
£ = s and n = 1; it falls from ~10% down to zero for N changing from 2 up to oo when £ = F 
{n = 1). 

In the latter case, ^ = F, for A'" = oo the weight of the first peak is 7.4% less than the exact 
value irh/a, and for high levels the weight tends to the exact value 7r7i|V'2fc(0)|^ = 7rh/a. 

Some deterioration of the results for AT = oo in comparison with the case N = 2 can be 
probably explained as follows. 

For N = 2 the straight line (/*(t) = 0, which serves as a leading term of expansion (5), is the 
unique classical path going from Xa = to = in time T, if T / 2TTk/uj {k = 0, 1, 2 . . .), 
that provides for the minimum of the action functional. When N 2 {in particular for N = oo) 
for any (even infinitesimal) T > there exists infinite (countable) set of classical paths (without 
a reflection, with one reflection, with two reflections, etc.) [24]. The influence of additional 
classical paths seems to grow for A^ — > oo. 
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3.24 




Figure 3: Anharmonic oscillators (24). Arrows point to exact values of even levels jl5j.j22j. 
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Table 1: Peak medians ill coinparisoii with, the exact Vcilues i/cxact 

[15] normalized as in (25) 



N 


2 


4 


10 


50 


oo 


Uexact 


3.1416 


2.9663 


2.9899 


3.2431 


3.4894 


yis 


3.08 


2.90 


2.82 


2.82 


2.832 


yiF 


2.79 


2.75 


2.84 


3.02 


3.145 


y2s 


3.12 


2.94 


2.86 


2.85 





3.3 Influence of continuum spectrum 

Above are considered the examples of potentials for which there exists only a discrete spectrum of 
eigenvalues of energy. For continuum spectrum the relation (2) is not applicable. It is therefore 
seems interesting to consider a well of finite depth, since this problem has both discrete and 
continuum ranges of a particle spectrum. 
In a well with the potential 

^{q)=-cosh-\q/j), (28) 

for distance unit a = Tr~^h{mUo)~^^^ in Eq.(8), already the first order of approximation using 
sinusoidal paths provides a reasonable accuracy. Plots of R\s{E') for various values of 7 are 
presented in Fig. 4. One can observe the growth of the number of bound levels E'^ < with 
growing depth of the well (growing parameter 7) , 
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E' 

Figure 4a: 
9 
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Figure 4d. 

It is important that not only approximations for ground levels are in good agreement with 
exact values [since for small q the potential (28) tends to a harmonic oscillator potential], but 
there is also a correct description of intermediate levels. The uppermost levels feel the influence 
of continuum. We get from (18) for n = 1 an asymptote for E' —0: 

Ru{E') ~ 7-i(-E')-^/' sin[77r-\-£;)-V2], (29) 

i.e. there appear oscillations with growing frequency and amplitude in the function R{E') for 
E' —>■ —0. The square root singularity R{E) ~ E"^/^ is typical for continuum spectrum, one 
can easily check this taking the integral in the left hand side of (2) for the evolution operator 
kernel of a free particle. In practice this singularity can be easily removed putting the system 
under consideration into a wide potential well, which transforms the continuum into a discrete 
spectrum. However, the distortion of uppermost levels is inevitable in this case also. 

3.4 A case of higher dimensions 

For oscillators with a potential 

if the dimension of the problem is i', we obtain from (18a) the following relation: 

Ris^E') = j{-isE')^''-^y^V2E'{l - x?)('^-i)/V,_i {js^E'xi{l - xlf'^) dxi 

D 
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where = Si^/r^ ^ and 5;^ is the area of the z/-dimensional sphere. 

Using the tables of integrals [25], we have from (30), leaving aside a constant factor: 

Ru2{E') <x E' Jl,^{i)- 



Ru4{E') oc {E'f{3Ji/^{0 - Jl/2(0^5/2(0}- 

Here ^ = V2jsE'; Uq = huj/2Tr. 

For multi-dimensional wells with perfectly reflecting walls with the potentials like 

4>{q) =0 for \q\<a 



(31) 



and 

we have 



0(9) = +00 for \q\ > a , 
RifuO^y"''^ j Ju-i{x)dx, 



where y = 2y^{ma'^E)^/^ /h. 

For u = 2 and = 3 it follows from (32) that 



and 



i?i/2 ocl - Jo(2V2y/7r) ; 

2\/2y7r 

i?i/3 OC y / Jo{x)dx - 2 Ji (2V2y/7r) 



(32) 



(33) 



The formulae (31) and (33) imply that medians of the peaks for ground levels are close 
to exact values of energy. The energy dependences of the peak weights obtained with these 
formulae also reproduce correctly the behavior of the squared module of the wavefunction: we 
have for oscillators |'02fc(O)P oc E^^'^^^'^ , and for wells, |V'2ifc(0)P oc E^^^^^"^ ■ 



4 Conclusions 

The numerical experiments and analytical examples considered in this work allow us to draw 
some conclusions. 

1. Replacing a path integral by its approximation with an integral of finite dimensionality 
may be an effective means for computing spectra of multi-dimensional quantum systems. 

It is remarkable that qualitatively, and even quantitatively (with accuracy ~ 10%), cor- 
rect information on the energy levels and on values of can be obtained already with 
approximations of low order, n = 1, 2, that is with a small amount of computing. 

2. Our preliminary computations, using Korobov numerical quadrature for n up to n = 16, 
show that the accuracy of method grows not faster than n~^/^. 
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3. It would be interesting to generalize the method for the case of many-body systems with 
account of particle statistics, probably with the help of Monte-Carlo methods, and also to 
apply this technique for computations of thermodynamic functions of non-ideal systems. 
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Addendum 

The text above is almost a literal translation of a Russian ITEP preprint 111-82. Only a few misprints 
arc corrected and a well known formula for £"„ in modified Poschl- Teller potential [a] is added after (28). 

During the two decades, that passed after the preprint was circulated, there appeared many papers on 
the subject of numerical applications of path integrals. Yet it seems that our results are still interesting 
because the technique of numerical analysis that we have undertaken, namely, using our formulae (18) 
and (19), was not explored by other workers for obtaining information on excited levels. Moreover, we 
found that the lowest orders in (18) and (19) are sufficient for obtaining reasonable results. Thus, this 
direction of research is still promising. 

Some historical references, missing in our text, are added below [a-e], as well as references to relevant 
papers [f-v] published after this report was issued in 1982. 

We thank Andrei Smilga for discussions that stimulated the translation of this paper into English 
and its posting in arXiv. S.B. is very grateful to H.Takabe for his kind hospitality at ILE, Osaka, where 
a major part of this version of the paper has been prepared, and to V.Zhakhovskii for useful hints on 
recent literature. 
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